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Qh' Abstract. Modern cosmology suggests that the Universe contains two dark components - 

dark matter and dark energy - both unkown in laboratory physics and both lacking direct 
evidence. Alternatively, a unified dark sector, described by a single fluid, has been proposed. 
Dissipation is a common phenomenon in nature and it thus seems natural to consider models 
dominated by a viscous dark fluid. We focus on the study of bulk viscosity, as isotropy 

f»o ', and homogeneity at large scales implies the suppression of shear viscosity, heat flow and 

diffusion. The generic ansatz £ oc p u for the coefficient of bulk viscosity (p denotes the 
mass/energy density), which for v = —1/2 mimics the ACDM background evolution, offers 
excellent fits to supernova and H(z) data. We show that viscous dark fluids suffer from large 
contributions to the integrated Sachs- Wolfe effect (generalising a previous study by Li & 

I"-- Barrow) and a suppression of structure growth at small-scales (as seen from a generalized 

Meszaros equation). Based on recent observations, we conclude that viscous dark fluid models 
(with £ oc p u and neglecting baryons) are strongly challenged. 
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1 Introduction 

The cosmological concordance model states that the Universe is spatially flat and approxi- 
mately 95% of its energy content is made up of an unknown dark sector. The remaining 5% is 
known: baryonic matter, electrons, photons and neutrinos. In the context of the concordance 
model a quarter of the dark sector behaves like cold dark matter (CDM), a pressureless com- 
ponent that clusters. The remaining dark stuff is called dark energy and is hold responsable 
for the current accelerated expansion of the Universe. 

Dark energy is commonly modeled either by a scalar field or a dissipationless fluid. The 
isotropy and homogeneity of the Universe at large scales, suggests that such a fluid description 
makes sense, at least as an effective description. Then dark energy has negative pressure, 
such that its equation of state p ~ — p today. Alternatively, one might view dark matter and 
dark energy as different manifestations of one single substance - unified dark matter. Its 
equation of state function w = p/p must be time-dependent in order to interpolates from the 
matter dominated epoch to the current accelerated expansion. 

In this work we study a single-fluid descriptions of the dark sector. As long as only 
the homogeneous and isotropic evolution (the background) is concerned, this class of models 
is indistinguishable from the concordance model. However, when it comes to perturbations 
differences become apparent at late times. We investigate in detail their contributions to 
the integrated Sachs- Wolfe effect, which probes most efficiently the large scales at late times, 
and the small-scale matter power spectrum, which probes the other side of the structure 
formation. We demonstrate that a generic single-fluid description of the Universe seems to 
be strongly challenged from both sides. 

A prominent candidate for this scenario is the Chaplygin gas, where p c = —A/p c , and 
its generalized version p gc = —A/p" [1]. For the case a = 1, this exotic equation of state is 
motivated by string theory, where the Chaplygin gas is interpreted as an effective description 



of a gas of D-branes in a D+2-dimensional space-time [2]. For the background evolution, 
a = corresponds to the A cold dark matter model (ACDM). 

Another popular candidate is a dissipative fluid with intrinsic bulk viscosity [3, 4]. Any 
real fluid shows dissipative phenomena and thus it is well motivated to include this aspect in 
cosmology as well. Typically shear viscosity is more important than bulk viscosity, however 
isotropy and homogeneity of the Universe at large scales, does not allow for shear. From the 
same argument diffusion and heat conduction cannot play an important role. Thus at large 
scales, bulk viscosity must be the dominant dissipative effect. For an expanding universe it 
gives rise to a negative contribution to pressure, p v = — £0, where is the volume expansion 
rate of the fluid and £ > is the coefficient of bulk viscosity. In order to study the cosmic 
expansion, one needs to specify £ = £(£). A common ansatz, which we will adopt below, is 

£ocp"[5]. 

The generalized Chaplygin gas (GCG) and the viscous dark fluid (VDF) have the same 

background dynamics in the one-fluid approximation [10]. However, the difference between 

these fluids appears at perturbative level. While the GCG is viewed as a dissipationless 

fluid, it has isentropic (the same entropy everywhere) perturbations 5p gc = (pg C /p gc )5p gc , 

the perturbative dynamics of the VDF is, by definition, nonadiabatic and it could also be 

non- isentropic, 5p v ^ (p v /p v )5p v . 

As GCG and VDF are equivalent w.r.t. their background evolution, both models can 
fit probes that are not sensitive to structure formation itself, e.g. SNIa data [6]. At the 
perturbative level one obtains, from the matter power spectrum and from the cosmic mi- 
crowave background (CMB) spectrum, very different predictions. The confrontation of the 
GCG model with the matter power spectrum data, provided by the Sloan Digital Sky Survey 
(SDSS) and the Two Degree Field Galaxy Redshift Survey (2dFGRS) data sets, discards 
a one-fluid GCG universe [7] due to strong oscillations in its theoretical power spectrum. 
However, it has been demonstrated that this problem is solved when a baryonic component 
is taken into account [8] or, if ad hoc entropy perturbations are included in the GCG pertur- 
bative dynamics [9]. While the former is a trivial solution, as the observed power spectrum 
corresponds to the visible matter and not the dark one, the latter can be seen as a motiva- 
tion to explore viscous (intrinsically nonadiabatic) cosmologies. On the other hand, the bulk 
viscous model does not show the same pathologies as the GCG, thanks to its nonadiabatic 
behavior [10, 11]. The GCG with a ~ agrees with CMB data [12], while, apparently, the 
VDF does not, due to a huge amplification of the integrated Sachs- Wolfe (ISW) signal [13]. 

The authors of [13] showed that, the evolution of the gravitational potential in a VDF 
model differs from the ACDM model at late times, implying a huge ISW effect. This is also 
found for the GCG [12] and for general unified dark matter cosmologies relying on a single 
scalar field [14]. However, these studies have been limited to fixed values of the cosmological 
parameters and it remains unclear, if the huge ISW effect could be avoided in a different region 
of parameter space. We study the dependence of the ISW effect on the model parameters of 
VDF models and GCG models. 

We also study the behaviour of sub horizon pertubations in the VDF and GCG models 
during the matter dominated epoch. We show that structure formation can be drastically 
affected in such cosmologies by comparing the growth of the unified dark matter perturbations 
with a typical CDM scenario. In other words, we investigate whether dark halos, the hosts 
of galaxies can form at all. 

In the next section, we compare the background evolution of the GCG and VDF models 
with ACDM. Section 3 is devoted to the study of linear perturbations. We provide an 



evolution equation for the study of the ISW effect and we obtain Meszaros-like equations 
for the evolution of sub-horizon perturbations in the VDF and GCG models. In section 4 
we derive quantitative results for unified dark matter cosmologies and conclude with some 
remarks and open issues in the final section. 

2 Background evolution 

In this work, we assume a spatially flat one-fluid description of the matter content of the 
Universe. This ansatz is expected to be appropriate at late times (thus radiation is negligible). 
At small scales we also neglect the effects of baryonic matter, which limits the precision of 
our discussion to the 10% to 20% level at small scales. 

The description of relativistic viscous fluids allows for a freedom in the choice of the 
comoving frame. Comoving observers could be comoving with energy transport (Landau- 
frame) or with particle number transport (Eckart-frame) . Both approaches are equivalent, 
but one has to make a choice. Here we adopt the Eckart formalism [15]. Then, the VDF bulk 
pressure is given by p v = — £0. Due to the second law of thermodynamics the coefficient of 
bulk viscosity £ > 0. The volume expansion rate O = u^ (Greek indices run from to 3, ";" 
denotes a covariant derivative) is obtained from the fluid velocity u A \ In a homogeneous and 
isotropic Universe, = 3H, where H is the Hubble expansion rate. With the ansatz 
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and assuming that the kinetic pressure p 
becomes [by means of H = Hq(p/pq) 1 ' 2 ] 
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The GCG model has a similar equation of state, p gc = —Ap$ (po/p) a , with A and a 
being dimensionless parameters. 

In a perfectly homogeneous and isotropic Universe, the GCG model and the VDF are 
equivalent, which is easily verified by the replacements a = —(y + ^) and A = 3i?o£o//?o- 
Hence, both fluids show the same time evolution. Instead of £o or A, it also is convenient to 
use the deceleration parameter q^. This correspondence can be established by 
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Once £o > 0, then qo < 1/2. The background evolution of the VDF and the GCG is governed 
by (a denotes the scale factor and ao = 1) 
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respectively. The existence of an early matter dominated epoch, H(a « 1) ~ a~ 3 ' 2 , is 

guaranteed for v < 1/2 and £o < Po/(3-Ho) for the VDF model and for a > — 1 and A < 1 in 

the GCG case. In order to obtain an accelerated epoch at late times (qo < 0), the parameters 

must obey £o > Po/(9Ho) an( A A > 1/3, respectively. The early and late time limits of both 

models are equivalent to the ACDM model. The only difference is the transition from the 

matter dominated phase to the accelerated epoch, which is given by the equation of state 

functions 
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The expressions in (2.5) are analogue to the ACDM one, 



§) =^f + i-n m0 , (2.8) 
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if we adopt q = ^f± - 1 (A = 1 - Q m0 ) and v = -1/2 (a = 0) for the VDF (GCG) model. 
These relations will be usefull in the next section in order to compare the perturbative 
dynamics of these models. 

3 Density perturbations of a dissipative fluid 

In this section we study the perturbative dynamics for the VDF and the GCG models. The 
differences between both models for an inhomogeneous Universe can be traced back to an 
inherent nonadiabatic behavior of the viscous model. In a sense, the VDF model can be seen 
as a nonadiabatic version of the GCG model. 

Let us start by considering the most general dissipative fluid with energy momentum 
tensor T^, including a dissipative contribution which is denoted by AT^. In the Eckart 
frame, the most general dissipative tensor is 

ATt = -£AT 6 ^ - V AT? V - kAT% v , (3.1) 

where £, ry and k are the coefficients of bulk viscosity, shear viscosity and heat conduction. 
For the homogeneous and isotropic background, only the bulk viscosity contributes to the 
cosmic dynamics. At first order, the heat conduction contributes only to the non-diagonal 
elements of AT^, and thus producing negligible contributions on superhorizon scales. The 
same happens with shear viscosity. In contrast to bulk viscosity, shear viscosity and heat 
conduction, influence the evolution of cosmological perturbations via spatial gradients. 

In the following we neglect heat conduction and shear viscosity also at the level of 
perturbations and thus the cosmic fluid is described by the energy-momentum tensor 

T% = pu»u v + pK + AT£ = pu»u v + pK - iv?. n h%, (3.2) 

where W 1V = g^ u + u^u u . More explicity, the background components of (3.2) are 

T° = -p, r° = T* = 0, Tj=pa8 i j =(p-^\s i j , (3.3) 



where H = — and the symbol (') means derivative wrt the conformal time r]. Latin indices 
run from 1 to 3. The effective pressure p e ff is the sum of an adiabatic component and the bulk 
viscous pressure (nonadiabatic). The VDF model is specified by p = and a dissipationless 
fluid is recovered with £ = 0. 

In the conformal Newtonian gauge the line element for scalar perturbations of an 
isotropic and homogeneous, spatially flat universe is 



ds 2 = a 2 (rj) [- (1 + 2(f)) dif + (1 - 2i/>) <%cfa;W] . 



The linear perturbations of the fluid 4- velocity are given by 
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For the linear perturbations of (3.2) we define the velocity scalar v, which is associated 
with the peculiar velocity by 5u l i = —kv/a, where k is the comoving wavenumber. The 
perturbed components of (3.2) read 
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5£ denotes the perturbation of the coefficient of bulk viscosity. The adiabatic speed of sound 
c| = (dp/dp)s- For dissipationless fluids, c| = p'/p' for the purposes of linear perturbation 
theory. For dissipative fluids in linear perturbation theory c| = {p' / p')^ = Q. 

As we neglect anisotropic stresses in our model, the spatial off-diagonal Einstein equa- 
tion implies (j) = ip. At first order, the (0-0), (Q-i) and the (i-i) components of the perturbed 
Einstein equation read (A = Sp/p) 



k 2 ^ - m^' - m 2 ^ = ^n 2 a, 
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The pressure perturbation 5p = cf^Sp + t5S, where 5S denotes entropy perturbations 
and r = (dp/dS) p . Below we assume that pressure perturbations do not give rise to spatial 
fluctuations of the entropy to baryon ratio. 

3.1 The integrated Sachs- Wolfe effect 

The ISW effect is a net change in the energy of a CMB photon as it passes through evolving 
gravitational potential wells. It can be computed by 

= 2 
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The integration is along the photon trajectory (h) from rj x (conformal time at recombination) 
to 770 (conformal time today). 

Combining the equations (3.9) - (3.11) into a single expression for the gravitational 
potential, we end up with 
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where S = <5£/£ can be considered as the relative perturbation of the coefficient of bulk 
viscosity. 

If we neglect the VDF contribution to the energy-momentum tensor, the right hand 
side of equation (3.13) vanishes and hence the resulting equation is the full evolution for the 
gravitational potential of an adiabatic fluid with an equation of state parameter w = p/p. 
The right hand side of (3.13) represents the influence of nonadiabaticity on ip. For the VDF 



model, we set cr s 



w 



and use the appropriate functions w v and H v . For the last term 



we need to know the functional form of £. If £ = £o(p/po) u its perturbation 5£ 
be related to the potential tp using equation (3.9). 
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3.2 Evolution of sub-horizon perturbations 

In the radiation era pressure suppresses the growth of structures. However, cold dark mat- 
ter, once kinetically decoupled from the plasma, starts to grow logarithmically on scales 
smaller than the Hubble horizon even during this epoch. Once the Universe becomes matter 
dominated (z C q ~ 3000) CDM can grow linearly in the scale factor. This scenario is called 
hierarchical structure formation as smallest structures form first and later on merge and grow 
to evolve into larger structures. 

The unified dark matter models studied in this paper have a matter-like behavior in 
the past, but do not necessarily provide a successful structure formation scenario. In order 
to study scales which entered the horizon sufficiently long before matter-radiation equality, 
we make use of the covariant conservation of the energy-momentum tensor (T^ M = 0). The 
first-order continuity equation reads 

A' - 3%A (w - 4 + i»v) ~ (1 + w + 2w v ) (kv + 3^') - mw^ip - S) = 0, (3.14) 



and the Euler equation is 



v' + 



n(i- 

w v k 



34 



3w v ) + - 



u\, 



IV- 



.A- 2 



+ w + w v 3% (1 + w + w v ) 



,, k(l + w) 
-tp + tt —ip + 



w v k „ 

£ + 



kc s 



H (1 + w + w v ) 1 + w + w v 1 + w + w v 1 + w + w v 



-A = 0. 



(3.15) 



For the adiabatic case there are many studies about the evolution of sub- horizon scales, even 
considering the possibility of energy other than matter or radiation [16] or modified theories 
of gravity [17]. However, the clustering properties of nonadiabatic CDM have not yet been 
considered in much detail. 



For the VDF model (w = c| = 0), we can simplify equations (3.14) and (3.15) and take 



the subhorizon limit of the Poisson equation to obtain 
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It is covenient to combine these equations to a single second-order differential equation 
for A and to use the scale factor a instead of conformal time. Hence, we obtain a Meszaros-like 
equation: 
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A = P(a), (3.19) 
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The function P(a) contains all contributions from the perturbation of the coefficient of 
bulk viscosity <5£. In the limit w v = we obtain the standard equation for CDM perturbations 
with the solution A cc j m oc a. 

The above equations are solved numerically. However, in order to obtain some analytic 
predictions for A, note that in the sub-horizon limit fc>H we find 
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If we also send w v to —1, the k 2 terms dominate and the equation is dominated by the first 
derivative term, thus one can expect exponential damping. 



4 Observational Constraints 

4.1 Supernova and H(z) data 

We employ a statistical analysis using recent H{z) [18] and the SN la constitution [19] data 
sets, in order to constrain the parameters of the background model. 

The confidence contours for a set of parameters {p} are obtained from the probability 
distribution function (PDF) 

P (p) = Be 2 , 
where B is a normalization constant. For a given sample, lets say SN, x 2 is defined by 

xi,(P) = E [ *' h(p) -/' (P)f («) 

i ' 

The quantities /if 1 and jj,° bs are the theoretical and the observed values, of the distance moduli 
and <Tj denotes their error for each data point i. For the H(z) sample we replace fj, by H. 
Hence, for the joint analysis we use \ 2 = x\n + Xh- 

Observational constraints on qo and v are shown in Figure (1). It displays the 2<r and 
3<7 confidence levels with best fit at (qo,v) = (—0.95,-3.2) with Xvi SC ous = 472.5. The 
dashed-red lines are age constraints for which the Universe is 13Gyr and 15Gyr old. The 
parameters for which the transition to the accelerated epoch occurs at Ztr = 1 and Zt T = 0.5 
are shown in the thin lines. We remark that these background results can be translated to 
the GCG model using the correspondences stablished in section 2. For the ACDM model 
(the horizontal line corresponding to v = —0.5) the best fit occurs at qo = —0.57 (vertical 
line) that means Q\ = 0.71. We obtain Xacdm = 472.9. Thus the latter x 2 is greater than 
for the viscous model. This occurs since the viscous model has an extra parameter. A model 
comparison by means of the Akaike information criterion, AIC = x 2 + 2fc with k being the 
number of free parameters [20], it becomes clear that both models are competitive with the 
ACDM model being slightly favoured (|AAIC| = 1.6). 

4.2 The integrated Sach- Wolfe effect from unified dark matter 

The CMB spectrum of anisotropies has been a key test for dark energy candidates as well 
as for modified gravity theories. It has been observed that UDM models suffer from an 
amplification of the ISW signal [13]. In general, for the GCG, unless a = 0, the acoustic 
peak to Sachs- Wolfe plateau ratio decreases for increasing a > 0. A similar conclusion was 
obtained for the bulk viscous fluid in [13]. However, the dependence of these results on the 
free parameters of the UDM models is still not clear and we adress this question now. 

We define a "quality" variable Q m to measure the difference between the ISW signal 
for some model m and the ACDM model, 

(AT\ m 

where positive (negative) values of Q stand for an enhanced (a reduced) ISW effect for the 
model m to ACDM. The signal (— ) can be obtained from (3.12), once we have calculated 
the gravitational potential ^> m from (3.13). A similar definition of Q was considered in [21]. 
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Figure 1. Observational constraints from SNla and H(z) data on the parameters of the VDF model 
(go and v). The fat dot indicates the best-fit model. Short-dashed lines denote 2a and 3er contours. 
Long-dashed (red) lines indicate age constraints of 13 Gy and 15 Gy, respectively. Thin lines denote 
the redshift of the onset of accelerated cosmic expansion. The cosmic expansion history of a ACDM 
model is obtained for v = —1/2 (horizontal line) with its best fit at qo = —0.57 (vertical line). 



The relevant modes for the ISW effect correspond to scales k < 0.003(/i/Mpc), that is 
the approximate scale where the Sachs- Wolfe C/ plateau begins in the CMB temperature 
anisotropy angular power spectrum. We shall plot the contours Q = 120%, 80%, 40% and 0% 
in parameter space and compare them to the background constraints obtained above. With 
this strategy we verify whether it is possible to conciliate the ISW effect contours close to 
Q = 0% (i.e. close to the ACDM model) with the "allowed" background model parameters. 

In order to estimate Q, we adopt a fiducial spatially flat ACDM model with parameters 
H = 72 km/s/Mpc and O m0 = 0.266, as suggested by WMAP-7. 

One can ask if the current measurements of the ISW effect are able to discriminate 
between different models. In other words, are the values Q = 40% and 80% or even Q = 120% 
acceptable? To answer this question we consider current and future estimations of the error 
bars of the CMB temperature and galaxy cross-correlation function C 9 t, that is currently 
used to measure the ISW effect [22] . Since the ISW effect is hard to measure one can currently 
discard only the models with Q > 100% (at 95%C.L.), corresponding to 25C g T/C g T > 1. 
Radio surveys in the near future will reduce the error up to a factor of 5 and thus should be 
able to improve the limit to the Q = 20% level (see figure 9 of [23]). 

4.2.1 A model for £(p) and its adiabatic counterpart 

With the ansatz £(p) = Co(p/Po) u the quantity H becomes 
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The evolution equation for the gravitational potential is obtained by combining (4.3) 
and (3.13). For the bulk viscous model (c 2 = 0) we solve it numerically and calculate Q v 



(see (4.2)) for various choices of the background parameters. The results are shown in figure 
(2). Q = 0% and even Q = 40% are in stark disagreement with the background contours 
(short-dashed lines) which are compatible with the constraints obtained in [24]. In additition, 
long-dashed lines display the age of the universe with 11, 13 and 15 Gyrs. The best-fit model, 
symbol • in Figure (2), corresponds to Q v = 120%. 




Figure 2. Additional CMB temperature fluctuations from the integrated Sachs- Wolfe effect are 
estimated by Q, see (4.2). Solid lines represent contours of constant Q v are shown in together with 
the 2t7 and 3cr contours and the age constraints of figure 1. VDF models with acceptable expansion 
history lead to at least a doubling of the ISW contributionwith respect to the WMAP 7yr best-fit 
ACDM model. 



For the GCG model we compute equation (3.13) with u> v = 0. Also we write Hasa 
function of A and a and for the adiabatic speed of sound we find 
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We observe, see figure (3), a small improvement, as the best fit model is close to the Q gc = 
80% line. However, both cases are discarded by the CMB data and this result agrees with 
[12-14]. 

In section 3 the perturbation of the bulk viscous coefficient (<5£) has been considered as 
a free function with its effects gathered by 3. Of course, a full perturbative analysis of the 
bulk viscous fluid should include this term. Let us for a moment take the freedom to treat 
3 as a free function of time neglecting the form imposed by (4.3). For the case 3 = cte = 
we observe that it is possible to conciliate the background constraints with the non-amplified 
ISW effect line (Q v = 0) as shown in figure 4. We remark that the background dynamics is 
exactly the same as before and the nonadiabatic contributions, except for the term 3, are 
still active on the r.h.s. of (3.13). The analysis shown in figure 4 reveals that 5 ^ is the 
source of the amplified ISW effect which has plagued the VDF model. For £ oc p u , we cannot 
regard 3 = as a solution to the problem, since 3 = occurs only if v = or A = 0. We 
know that there are density fluctuations in the Universe, thus A/0. On the other hand, 
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Figure 3. As figure 3, but now for the GCG model. The free parameters of the GCG model are 
A and a, as defined in the text. Dashed lines represent the 2a and 3a contours of the fit to SNla 
and H(z) data, with best-fit denoted by the fat dot. From left to right, the solid lines are contours 
of constant ISW contribution Q gc = 0%, +40%, +80% and +120%. For the GCG the ISW effect is 
slightly less pronounced compared to the VDF models. 



assuming v = leads to very different background dynamics, which we study below as a 
particular configuration of the VDF model. 




-0.975 -0.95 -0.925 -0 . 

q0 
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Figure 4. As figure 2, but now the perturbation of the coefficient of bulk viscosity is arbitrarily set 
to zero, St; = 0. ^From top to botton the solid lines are the corresponding Q v = +120%, +80%, +40% 
and 0% contours. Long-dashed (red) lines display the age of the universe with 11, 13 and 15 Gy. 
Now, VDF models fit the background and do not show an enhanced ISW effect. However, there is no 
physical motivation to put <5£ = in the context of VDF models, unless v = 0. 



We have also verified that extremely large negative (positive) values for v(a) do not 



If 



produce a large amplification in the ISW effect. Concerning the GCG, this limit of the 
parameter a had been found before in [25] but the correspondence with v had not yet been 
established. Large values for the parameter v also agree with the analysis using the matter 
power spectrum [10]. This range for the parameter v{pi) implies a step-transition of the 
background evolution from a CDM phase to a deSitter one as discussed in [26]. 

4.2.2 A constant coefficient of bulk viscosity 

The previous considerations suggest to study the case of a constant coefficient of bulk viscosity 
[y = 0) in more detail. Now, the VDF has only qo as free parameter. At perturbative level 
there are no contributions from H but the rhs of (3.13) is non- vanishing. Figure 5 shows the 
PDF for go with the values Q v = 120%, 80%, 40%, 0% and constraints from the age of the 
universe. The line Q v = 0% is within the 2a region, but leads to a Universe younger than 
13 Gyrs. In order to satistfy the age constraints Q > 30%, which can be tested in the near 
future [23]. 
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Figure 5. PDF for the case v = with best fit at go = —0.46. The short-dashed lines denote the 
la and 2a regions. The age constraints {13Gy and 15Gy) are shown by long-dashed lines. Solid lines 
represent, from left to right, Q v = 120%, 80%, 40% and 0%. The regions of a small enhancement of 
the ISW effect is in conflict with the age of the Universe. 



4.2.3 Mimicking the ACDM background evolution 

For v = -0.5(a = 0) the VDF (GCG) and the ACDM models have exactly the same 
background evolution. Hence, the ISW contribution from nonadiabatic perturbations can 
be quantified. Note that identical background solutions can be achieved by three different 
models: i)the ACDM scenario, ii) a two-fluid model consisting of pressureless, dissipationless 
matter and a dissipationless fluid with EoS p = —p, iii) the VDF(GCG) with v = — 0.5(q = 0) 
which is also equivalent to a fluid with a negative constant pressure. On the other hand, these 
models have distinct perturbative dynamics, namely: i) there are no perturbations from A, 
ii) could have nontrivial but adiabatic perturbations and iii) has nontrivial and non-adiabatic 
perturbations. The GCG is an example for ii). 

The result for the VDF with v = —0.5 is shown in the left panel in Figure 6. The 
nonadiabatic contributions of the VDF are responsible for putting the Q v = 0% line outside 
3<t CL. However, if we neglect the contribution from <5£, the Q v = line is within the la CL, 
right panel in Figure 6. 
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Concerning the GCG with a = 0, the PDF for A parameter is shown in Figure 7. This 
particular case behaves very similar to the ACDM and our result agrees with [12]. 
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Figure 6. PDFs for the VDF with v = —0.5. The best fit occurs at qO = —0.64. Left panel shows 
the result considering the full evolution while in the right panel the pertubation S was neglected. The 
age constraints (13 Gy and 15 Gy) are shown as long-dashed lines. Solid lines represent, from left to 
right, Q v = 120%, 80%, 40% and 0%. 
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Figure 7. PDF for the GCG with a = and best fit at A = 0.76. The age constraints (13 Gy 
and 15 Gy) are shown as long-dashed lines. Solid lines represent, from left(right) to right(lcft), 
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120%, 80%, 40% and 0%. This GCG model does not suffer from an ISW overproduction 



To summarize the observational constraints of sections 4.1 and 4.2, we have seen that 
generic VDF models that are excellent fits to SNIa and H(z) data sets, give rise to a large 
ISW contribution to the CMB temperature angular power spectrum and are thus excluded. 

4.3 Structure formation on small scales 



In the standard CDM structure formation scenario small-scale perturbations start to grow 
oc a when the universe becomes matter dominated, at a redshift z cq . Before z eq , even if 
the wavelength of the perturbation is larger than the Jeans length rapid expansion prevents 
the growth of structures. Hence, before we can study the process of strucutre formation for 
the VDF it is essential to establish the time at which the universe beconres VDF dominated 
and we thus include the radiation fluid in our analysis. With the inclusion of radiation 
the dynamics of the GCG remains the same. However, since the expansion rate becomes 
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H = [^p(p v + Pr)] 1 ' 2 the background dynamics of the VDF is severely changed at early 
times. The fractional density for the VDF is given by the numerical solution of 



a— h 3S^ 

da 



^+^¥ 



1/2 



0, 



(4.5) 



where £ = 9-ffo£oPc , p c is the critical density and O r o = 8.475 x 10 . The new model 
parameter £ is related to the deceleration parameter qo approximately by 



Qo 



5' 



PvO + PrO 



fr 



0- 



(4.6) 



The fiducial ACDM model adopted in section 3 has the matter-radiation equality occuring at 
z C q = 3137. For a VDF plus radiation the equality is a function of the model parameters and 
will be denoted by z* q . As shown in figure 8 for parameters values within 2a CL, z* q > z cq . 
Hence, sub-horizon VDF fluid perturbations start to grow before typical CDM perturbations. 
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Figure 8. Redshift of matter-radiation equality for the VDF (left) and GCG (right). Countours of 
constant z* (see text) are shown in the model parameter space of the VDF (left) and GCG (right). 
Dashed lines denote the constraints from SN and H(z) data as presented in figures 1 and 3, respectively. 



We solve equation (3.19) with intial conditions A v (z cq *) = 1 and -^(^eq*) = 1 arid 
compare with the CDM evolution A ct j m oc a with the same initial conditions, however cal- 
culated at z cq , A cdm (z eq *) 
function in (3.19) become 



1 and — ^r^^cq*) = 1- The Hubble rate and equation of state 



Hn 



Q v + Q, r0 a~ 



UK 



1-2^0 



(n v + Ooa- 4 ) 1 / 2 ^, 



(4.7) 



with Q v being determined from (4.5). 

We consider modes which give rise to cluster (subgalactic) size structures k ~ 0.2Mpc~ 
(k = 10 6 Mpc~ ). We assume that for these modes the nonadiabatic Meszaros equation is 
valid up to the onset of non-linear evolution at z n \ = 3(60 ± 20) [27]. Soon after z n \, a large 
fraction of the matter collapses into gravitationally bound objects. Nonlinear effects lead to 
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a further modification of the final (at z = 0) power spectrum. The study of them is beyond 
the scope of this work. Figure (9) shows the growth of perturbations, considering the best 
fit model obtained above, for k = 0.2 - 0.3Mpc _1 (k = 10 6 Mpc _1 ) in left (right) panel. Also 
the CDM growth is shown as short-dashed line. If we consider the full evolution of equation 
(3.19) including the term 3 (botton lines indicated by 5£ ^ 0) we observe a large growth 
suppression after a redshift z ~ 6(a ~ 0.14) for k ~ 0.2Mpc~ and z ~ 200(a ~ 0.005) 
for k = 10 6 Mpc~ . Indeed, the dominant contribution in the terms proportional to k 2 A 
and k 2 A' comes from 3 and, consequently, at late times the density contrast A will decay 
rapidly. On the other hand, similarly to the ISW effect results, the perturbative dynamics is 
well behaved if <5£ = (upper lines). 

The GCG perturbations do not suffer any kind of suppression and will behave exactly 
like the CDM ones since it obeys the standard adiabatic growth equation, w v = in (3.19), 
with solution A gc oc a. At the same time, for any GCG configuration the transition to the 
accelerated expansion phase occurs after z n i and the growth of perturbations is not suppressed 
by this effect. 



k=10 6 Mpc" 




<5f=0 



0.005 0.01 0.015 0.02 



Figure 9. Left panel: Growth of sub- horizon perturbations in CDM (short-dashed) in the ACDM 
model and of the VDF for k = 0.2Mpc~ (long-dashed) and k = 0.3Mpc~ (solid). The upper lines 
for the viscous fluid have 6£ = 0, while the bottom ones have 8£, = is£A. Right panel: The same for 
k = 10 6 Mpc~ . Generic VDF models supress structures on subgalactic scales exponentially. 



5 Conclusions 

The main idea behind the unification scenario is to reduce the dark sector to one component 
instead of dark energy and dark matter. This component should, at cosmological scales, 
reproduce both the structure formation process and the current accelerated expansion of the 
universe. The former condition seems to be the main challenge for such models. 

We have compared the ISW signal of UDM models with the ACDM prediction. Figures 
(2) and (3) are in agreement with previous results, where the background-prefered model 
parameters of the VDF and the GCG imply an unacceptably large amplification of the ISW 
effect. In fact, we confirm and quantify the findings of [13] for a wide range of parameters. 
Although models with u = cannot be ruled out by current data, they nevertheless show a 
significant amplification of the ISW effect that will be detectable in the future. 

This tight constraints can be seen as an evidence that either bulk viscous effects do 
not play a role in the cosmic dynamics or that the phenomenological ansatz £ ~ p v is not 
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appropriate. Since the intensive thermodynamic variables are functions of the extensive ones, 
a possible alternative is to describe the bulk viscous pressure in terms of energy density and 
entropy, i.e. p = p(p,S) and £ = £(p, S). This could imply a well behaved perturbative 
dynamics and alleviate the ISW problem of such fluids. Recently a microscopic model for 
the cosmic bulk viscosity has been introduced as a dark energy candidate in [28]. This "dark 
goo" model shows good results when compared with the matter and CMB power spectrum. 
It could be interesting to extend this model to the context of unified dark matter where an 
estimation of the mass of the dark particle can be achieved. An important lesson from "dark 
goo is that more realistic viscosity coefficients show a complicated dependence on the energy 
density and cannot be written as £ oc p v . 

On the other hand, we note that bulk viscous pressure represents a small negative 
correction to the positive equilibrium pressure. Here we have admited the viscous pressure 
to be the dominating part of the pressure. This is clearly beyond the established range of 
validity of conventional non-equilibrium thermodynamics and non-standard interactions are 
required to support such an approach [29]. Hence, viscous cosmologies based on the Israel- 
Stewart theories [30] can also be considered. Recently, a qualitative analysis of such causal 
transport theory has been performed in [31]. 

We have studied the evolution of sub-horizon scales during the matter dominated epoch. 
In the standard CDM model the linear growth of small-scale perturbations gives rise to dark 
halos hosting galaxies. Concerning the unified scenario, we find a modification of the redshift 
of the matter-radiation equality. As shown in figure 8 the prefered parameter values for the 
UDM models are compatible with z* q > z eq . Hence, UDM perturbations start to grow earlier 
than CDM perturbations. The GCG perturbations follow the CDM growth A gc oc a until 
z n i and consequently, only the amplitude of the perturbations will be different. On the other 
hand, the VDF perturbations grow in a different way following a nonadiabatic Meszaros- 
like equation derived in section 3. In general, the evolution of A is scale-dependent and 
deviates significantly from ACDM. The most important effect at late times is the dominance 
of nonadiabatic contributions causing A to decay rapidly. 

Despite the different evolution of viscous-matter perturbations and standard CDM per- 
turbations on subhorizon scales, the ad hoc assumption 5£ = can alleviate the A growth 
suppression. Since VDF perturbations start to grow before z cq , their amplitudes are of similar 
size as in the ACDM case. In other words, the growth before z n \ offsets the late time growth 
suppression. But this assumption is not an acceptable solution and unless other effects like 
shear viscosity or a very different ansatz for £ would lead to qualitatively different results, the 
VDF models are ruled out. This can be interpreted as a complementary probe that viscous 
cosmologies based on the Eckart formalism (with £ oc p u ) are strongly challenged as potential 
contenders for a general relativity-based description of the cosmic medium. 

The inclusion of a baryonic component in the system provides a more realistic model 
when compared with the one-fluid approximation adopted here. In this case, the background 
dynamics of the GCG fluid remains the same while the VDF will behave differently. How- 
ever, as baryons represent a small fraction of the cosmic energy budget we do not expect a 
significant influence on the evolution of the unified dark sector. 

To conclude, the dissipative UDM models considered in this work have severe problems 
to describe the observed cosmic structure on largest (enhanced ISW effect) and on smallest 
scales (overdamping due to dissipation). 
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